Ground state of graphite ribbons with zigzag edges 
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We study the interaction effects on the ground state of nanographite ribbons with zigzag edges. 
Within the mean-field approximation, we found that there are two possible phases: the singlet 
superconducting (SS) phase and the excitonic insulator (EI). The two phases are separated by a 
first-order transition point. After taking into account the low- lying fluctuations around the mean- 
field solutions, the SS phase becomes a spin liquid phase with one gapless charge mode. On the 
other hand, all excitations in the EI phase, especially the spin excitations, are gapped. 
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I. INTRODUCTION 

After the discovery of low-dimensional materials such 
as fullerenes and carbon nanotubes, the research on the 
sp 2 network systems has been attracting much attention. 
The nanographite ribbon is one of the most simple and 
fundamental fragments of the sp 2 network and represents 
a new class of mesoscopic systems. In this system, the 
boundary regions play an important role so that the edge 
effects may influence strongly the 7r-electron states near 
the Fermi surface. 

There are two basic shapes of regular graphite edges 
- zigzag and armchair edges (see Fig. |l|). The study 
of electronic states of hydrogen-terminated graphite rib- 
bons reveals that the ribbons with zigzag edges possess 
partly flat bands at the Fermi level, which correspond 
to the electronic states localized in the near vicinity of 
the edgesH'BS. Specially, the highest valence band and 
the lowest conduction band are always degenerate at 
kao = 7r with the lattice spacing ao ~ 2.46 A. (Hereafter, 
we will set ao = 1.) The localized edge states are of 
special interest because of their relatively large contribu- 
tion to the density of states (DOS) at the Fermi surface, 
which results in the CurieJike temperature dependence 
of the Pauli susceptibility^ and zero-conductance reso- 
nances in the nanographite ribbon junctionso. It was re- 
ported that the zigzag ribbons do not undergo bond alter- 
nations along the ribbon axis for reasonable strength of 
electron-phonon interactions, because of the non-bonding 
character of the edge statesO. In other words, the partly 
flat bands are stable against the Peierls instability. In 
addition, the flat edge states exist not only in the single- 
layered zigzag ribbons but also in the stacked layers of 
zigzag ribbons in a manner of the AB stackingB, in which 
half of the carbon atoms of one ribbon are located di- 
rectly above the center of each hexagon on the neighbor- 
ing ribbons. 

In this paper, we are interested in the low-energy 
physics of the graphite ribbons with zigzag edges. In 
that case, it suffices to consider the edge states only. A 
simple power-counting indicates that all four-fermion in- 
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FIG. 1: The structure of graphite ribbons with zigzag edges. 
The carbon atoms are located at the corners of each hexagons. 
The rectangle with the dashed line is the unit cell. 



teractions are relevant operators around the free-fermion 
fixed point (see section III). Therefore, we expect that the 
electron-electron interactions will change the free-fermion 
picture drastically. Based on the renormalization group 
(RG) analysis, we propose a model Hamiltonian with 
an 0(4) symmetry to describe the low energy physics 
of this system. A mean-field treatment of this Hamil- 
tonian results in two possible zero temperature phases 
- SS and EI phases separated by a first-order transition 
point. For the EI phaseQ, the singlet and triplet exci- 
tonic order parameters can coexist due to the underlying 
0(4) symmetry!! As a result, there is a degeneracy for 
the onset of the charge density wave (CDW) and spin 
density wave (SDW) ordering within the mean-field ap- 
proximation. However, in one dimension, the low-lying 
fluctuations around the mean-field solutions are so strong 
that the long range orders obtained by the mean-field 
theory are destroyed and the true ground states exhibit 
the algebraic (SS) or short-ranged (EI) orders. Conse- 
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FIG. 2: The schematic figure of the energy band dispersion 
near E = 0. A z is given by Eq. (^|). The band indices 1 and 
2 indicate the two lowest bands. 



quently, the SS phase becomes a spin liquid with one 
gapless charge mode, whereas the EI phase turns into an 
insulating phase in which all excitations are gapped and 
the broken 0(4) symmetry is restored. 

The rest of the paper is organized as follows: In section 
II, we give a description of the system and discuss the 
action describing the dynamics of the flat edge states. We 
present the renormalization group (RG) analysis of the 
action in section III and propose an effective Hamiltonian 
which describes the low energy physics. The mean- field 
theory and the derivation of low-energy effective actions 
are given in section IV and V, respectively. The last 
section is devoted to the discussions and conclusions. 



II. THE MODEL SYSTEM 

The structure of the graphite ribbon with zigzag edges 
and the schematic figure of its energy band dispersion 
near E = are shown in Figs. |l| and ||, respectively. The 
almost flat bands appear within the region < \k\ < n. 
For an H-terminated single-layered zigzag ribbon, the dis- 
persion relations of the two lowest bands close to fc = tt 
obtained ffcan the tight-binding model has the approxi- 
mate formtffl 
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where = 2cos(|), t is the hopping matrix clement 
and N is the number of zigzag lines (see Fig. [j]). Eq. 
(|I|) indicates that the two lowest bands are degenerate 
at the Fermi point k = fco — ir. Around that point, 
£i(2)k ~ ±2Nt\k — kol 1 ^^ 1 - The energv gap to the next 
band A 2 as shown in Fig. || is given bjoEl 

~(N- 1)tt 



A z = At cos 



(2) 
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Since the energy scale we are considering is much lower 
than A 2 , only these two lowest bands are involved and 



the electron operator can be expanded around the Fermi 
point fcg 



*a(£) 



ikox 



i/)i a (x)ui(y) + e a pi/>lp(x)u2(y) , (3) 



such that ip aa \0) = 0, where a = 1,2 is the band in- 
dex, and a = ±1 corresponds to spin up and down, re- 
spectively. (We assume that the x-axis is parallel to the 
ribbon axis.) The -0-fermions describe the low-energy de- 
grees of freedom. 1/1.2(2/) are re& l functions which satisfy 
the orthonormal condition: J dy u a (y)uj,(y) — S a b- Un- 
der the operation of reflection about the line along the 
ribbon axis at the middle of the ribbon, denoted by a 
unitary operator 1Z, we have IZui — ui. (Note that this 



results in u\ 



•) 



Up to the four-fermion interactions, the most general 
form of the action describing the dynamics of -0-fermions 
in the imaginary-time formulation is given by 
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drdx Li 



where 



h 
Li 



% / A A §^L(-^+^nvw, (4) 



y (pi + Pi) + (J2P1P2 + 4.g 3 Ji ■ J2 
^(NtNi + H.c.) . 



(5) 



Here ip aa (k) = J dx e tkx ip aa (x), m > 1 (m = N — 1 



for the single-layered graphite ribbons) , v 
UV cut-off for the momentum, and 
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2Nt, A is an 
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In Eq. (|), we have replaced fc — fco by fc. p a , J a , and 
N a defined above are, respectively, the charge density, 
the spin density, and the singlet Cooper pair in the band 
a. g\ is the intra-band interaction, and the inter-band 
interactions are described by the density-density inter- 
action g 2 , the exchange interaction g 3 , and the singlet 
Cooper-pair tunneling 94. (Note that the triplet Cooper- 
pair tunneling term vanishes in the present case because 
of the Fermi statistics.) We have to emphasize that, re- 
gardless of the graphite ribbons being H-terminated, not 
H-terminated, single-layered, or stacked-layered, Eqs. (Q) 
and (@) describe the dynamics of the flat edges states as 
long as they exist. The short distance structure only af- 
fects the values of the parameters gi (i — 1, • • • ,4), v, 
and m, which will ultimately determine in which phase 
the system is located. 

To obtain the values of the couplings gi (i = 1, ■ ■ ■ , 4), 
we consider the short-ranged electron-electron interac- 
tions described by the interacting Hamiltonian 
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Hint = ~ / d 2 x Y d 2 x 2 : tft(£i)*J(£ 2 )V(2i - x 2 )^0{x 2 )^a(xi) : . (7) 
By inserting Eq. (0) into Eq. (]?]), we obtain 

31 = ^(0) , g 2 = -Ui(0)+ l ??p , -2g 3 = U 2 (0)=g 4 , (8) 
where U l {q) = / dx C^V^x) (i = 1,2), and 

1^.(0:1-12) = y dy 1 dy 2 u\(yi)V(xi - x 2 )u\{y 2 ) , 

V 2 (xi -x 2 ) = / dyidy 2 u\ (yi)u 2 (yi)V(xi - x 2 )m {y 2 )u 2 (y 2 ) . 



Especially, for a Hubbard-like interaction V{x\ — x 2 ) = 
Vq5(x\~ x 2 ), we have gi = U = —2g 2 = —2g3 = .94 where 
U = Vo J dy uf(y). We see that g\, g^ > and g 2 , 33 < 
when the interactions between electrons are repulsive. 

The symmetries of the action / is the charge U(l) 
(denoted by U c (l)) where the V'-fermions transform as: 
"01a ~~ ¥ e %e ipi a and ip 2a — > e i/) 2a , the spin SU(2), and 
the Z2 (particle- hole) where the ^-fermions transform as: 
ipi a «-> 4> 2a . In this paper, we restrict ourselves to the 
undoped case. The chemical potential term breaks the 
Z2 symmetry and thus will not be generated by renor- 
malization. 



III. RENORMALIZATION GROUP ANALYSIS 

We first show that L\ is a relevant perturbation to the 
free-fermion action Iq. By integrating out the fast modes 
with A/s < \k\ < A, we find that Jo is a fixed point of 
the following RG transformation 

k — >fe/s, w -> ws" m , 4> aa , (9) 

with £ = s'™ -1 / 2 . Under the RG transformation (||), the 
couplings in L\ transform as g\ — s m ~ 1 gi (i = l,---,4). 
We see that in one dimension all couplings are relevant 
and have the same scaling dimension around the free- 
fermion fixed point Iq . To tell which terms dominate the 
low energy physics, a controllable approximation to or- 
ganize the quantum fluctuations is necessary. Here we 
adopt the e-expansiorJa. That is, we extend the spatial 
dimension from 1 to d and then use e = d u — d as the ex- 
pansion parameter to compute the RG functions, where 
d u = m is the upper-critical dimension. For small TV 
(and thus small m), the e-expansion may be reliable. For 
larger N, it is hoped that there is no qualitative differ- 
ence in physical properties from small N to large N as 
long as N is finite. 

Next we would like to examine the quantum fluctua- 
tions up to the one-loop order. Defining the dimension- 
less couplings: A; = A~ e cgi/v (i = 4) where c = 



Sn-i/[2(2tt) n ^ 1 ] and Sd is the area of a d-dimensional 
sphere, and calculating the one-loop corrections to the 
coupling constants g[ (i — 1, • • • ,4), we obtain the one- 
loop RG equations within the e-expansion: 

dXl - P \ X 2 X 2 
~df ~ 1 4 > 

dX2 - f \ X 2 IX 2 X 2 
— — — EA 2 — A 2 — 0A3 — A 4 , 

—zr = £ A 3 — 2A 3 (A2 — A3) , 

= eA 4 - 2A 4 (Ai + 2A 2 ) . (10) 

Eq. ([ll]) is solved numerically up to a scale £*, at which 
point the largest coupling X m ax — max{Ai(i*)} satisfies 
U -C Xmax <C 1, where U = A~ e cU/v. This allows us 
to ignore the higher-order terms [0(Af )] in the RG equa- 
tions. As U — > 0, I* — > 00, and we need only analyze 
the asymptotic large I behaviors of Eq. ([l0]). We find 
that for e > 0, A2 and A4 become divergent first while 
the value of Ai approaches zero and that of A3 remains 
small when A2.4 = 0(1). 

Based on this analysis, we propose that the low-energy 
physics of the graphite ribbons with zigzag edges is de- 
scribed by the following model Hamiltonian: 

/dk 
- e^Uk)^aa(k) (11) 

+ J dx [g(NiN 2 + H.c.) - g Pl p 2 ] , 

where <?,<? > and = u|/c| m . The symmetry of 
H is U c (l)xSU(2)xSU(2)xZ 2 where the SU(2)xSU(2) 
symmetry corresponds to the independent spin rota- 
tions in each band. The enlarged symmetry (SU(2)— > 
SU(2)xSU(2)) arises from the suppression of the 33 term 
under the RG transformation. The values of g and g de- 
pend on v and the initial values of <7j's (i = 1, • • • 4). In 
the following, we shall treat the coupling constants g and 
g as free parameters and study the possible phases of the 
Hamiltonian (|TT|) . 
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IV. MEAN-FIELD THEORY 



II 



The g term in the Hamiltonian (jll|) describes the pro- 
cess of singlet Cooper-pair tunneling between two bands 
and favors the singlet superconductor when it becomes 
divergent, whereas a strong g term enhances the fluctua- 
tions of electron-hole pairs or excitonic ordering between 
band 1 and band 2B. The latter can be represented by 
the singlet or triplet excitonic order parameters according 
to the angular momentum it carries. Motivated by this 
observation, we define the following order parameters: 

Osa(x) = gip a -\(x)lpal(x) , (12) 

With the help of the Hubbard-Stratonovich transforma- 
tion, the corresponding action of the Hamiltonian ( |ll| ) in 
the imaginary-time formulation can be written as 



S 



duj dk ~ + ~ 

Z7T Z7T 

■J drdx (C1+C2+C3) 



(13) 



where 

£3 



^OslOs2 - d s l1p2f1p2l - O s2 4>lf1pll + H.C. 

^|$ 5 | 2 - (iQeaptpiatop + H.C.) , 
9 v J 

2 / ~ 

-|*t| 2 - [3>\ ■ e aX il>2\{cr) a ^ip + H, 



The mean-field treatment of the action S ([T^) is to 
assume the presence of bosonic mean fields, neglect the 
fluctuations of order parameters, and finally, integrate 
out the fermionic degrees of freedom to derive the effec- 
tive potential. Now we consider the mean-field ansatz: 
O sa = (O sa ), $ s = (<£> s ), and * t = (* f ). The effective 
potential is given by 



V 



-(O sl O s2 + C.c.) + i(\^ 
9 9 
F(uj,k) 



l* t | 2 ) 



dio dk 
2tt 2tt 



In 



F (u,k) 



(14) 



where Fq(lj, k) = {oj 2 + e|) 2 , and 

F(«,fc) = (cj 2 + etf + \0* sl O* s2 + A T A\ 2 

+ (c 2 + e 2 fe )(|O sl | 2 + |O s2 | 2 + 2AtA) . 

Here A = ($ s , $t) T is the excitonic order parameter. In 
Eq. (|i"4|), we have set V = for free fermions. 

The ground state is determined by the minimum of the 
effective potential V and the solutions of the mean-field 
equations are given by 



A 



9/$ 



FIG. 3: The mean-field phase diagram. Region I is the EI 
phase and region II is the SS phase. The point A (g = g) is a 
first-order transition point 



where A s 
R ( 1 I 

D V2m' 2 



C m V m-lgm-l with C n 



and a ri 



/(4m7r), and 



A % 



O, 



= O s 



(16) 



where Ao is a real vector with the fixed length: A e = 

Here B(x,y) is the beta 
Ttih correspond to the 



Osi = 0% = A s e 



, 



(15) 



V A A = CmV —ig- 

function. The solutions ( |l5| ) and 
SS and EI phases, respectively. Especially, in the EI 
phase, 7^ can coexist with 3? t 7^ 0. It is well 
knowrollj that the singlet ($ s ) and triplet (4>t) exci- 
tonic orders are accompanied by the appearance of CDW 
and SDW. In our case, the onset of both ordering is de- 
generate on account of the 0(4) or SU(2)xSU(2) sym- 
metry of the action ([l3]). The weak interactions which 
are neglected in the Hamiltonian ( pi] ) will lift this de- 
generacy. A residual Coulomb interaction like the g^ 
term favors the triplet excitonic order, while the electroa,- 
phonon interactions favor the singlet excitonic ordered. 
Note that naively the hidden symmetry of the EI phase 
is U/(l)xO(4), where the U/(l) symmetry corresponds 

to the transformation: ^p aa — > e % ®il)g a . However, it is not 
a symmetry of the Hamiltonian ( |ll|) and we will see later 
that after taking into account the fluctuations the value 
of 9q cannot be arbitrary. 

The point g = g corresponds to the first-order phase 
transition because at that point the ground-state energies 
of the SS and EI phases are equal and the correspond- 
ing order parameters do not vanish. Thus, the values 
of the order parameters change discontinuously from the 
SS to EI phases or vice versa. This is consistent with the 
RG flow obtained from the e-expansion. Since the RG 
equation (|l(]) docs not exhibit any IR stable fixed point, 
the phase transition cannot be the second-order one. To 
sum up, the zero-temperature phase diagram within the 
mean-field approximation is shown in Fig. ^. 



V. LOW-ENERGY EFFECTIVE ACTION 

To examine the stability of the mean-field solutions ob- 
tained above, we have to consider the fluctuations around 
the mean-field ground state. 

In the SS phase, O s \ and S 2 can be parametrized as: 
6 si = A s e i01 and 6 s2 = A s e"^ 2 . On the other hand, 
the excitonic order parameter A becomes the ordinary 
bosonic field. Inserting these into the action ([l3]) and 
integrating out the fermionic fields and A, we obtain the 
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low-energy effective theory in the SS phase 
1 



Ko 
2 

4 



-(dr^+y +v s (d x( f> + y 



(17) 



i 



-(d T (j>-Y + v s {d x( j>-Y 



Ho cos 4>- 



where d>+ = d>i ± ( 



fa, K = K Q = ^j^-Vambm, v s = v s = 

l— — i 
mW^-u^As m , and b rn oc (A 2 /A s ) 1_ ~ 3> 1 is a non- 
universal constant. The (/> + -ficld carries charge two and 
spin zero while the quantum numbers of <fi- are Q = 
and 5 = 0. The cosine term in C s opens a gap in the <p- 
sector because Kq > (87r) _1 . The charge U(l) symmetry 
forbids the terms like cos ((3(j)+) or sin (/3(f> + ). Therefore, 
the <p+ sector remains gapless. After taking into account 
the fluctuations, the mean- field gap A s becomes the spin 
gap and the SS phase is, in fact, a spin liquid (SL) phase 
with one gapless charge mode, which is sipilar to the 
C1S0 phase in the two-leg Hubbard laddeiO. Besides, 
the long range SS order turns into the algebraic one due 
to the gapless 0+-field and the SS fluctuations are en- 
hanced compared with the free fermions. 

In the EI phase, the excitonic order parameter A can 
be written as: A = A e $e ie where $ is a real vector 
satisfying $ T <I> = 1 and transforms with the fundamental 
representation of 0(4), whereas O a \ and O s i are ordinary 
bosonic fields. The quantum numbers of the #-field are 
Q = and 5 = and the ones carried by the $ sector 
consist of (Q, S) — (0, 0), (0, 1). The low-energy effective 
theory in this case can be obtained by inserting the above 
parametrization into the action ([l3]) and integrating out 
the fermionic fields, O s i, and S 2, and then we have 



S 



EI 



J drdx (Co + £$) , 



where 
Co 
C$ 



K 
~2 

Po 

2 



-(d T e) 2 + Ve {d x ef 



(d T ^f+V (d x $f 



+ A cos (29) , (18) 
(19) 



Here K = p = 2(m - l)\/a m 6 m , v e = v = 
ii^/^umAe m , and b m oc (A z /A e ) 1_ m ^ 1 is a non- 



universal constant. In Eq. (18), we keep the leading 
cosine term only. Because K >(2it)~ 1 , the cosine term 
in Ce is a relevant operator and the 9 sector acquires 
a gap. Beisdes, (9) is pinned at some value 9$ which 
depends on the short distance physics. The dynamics of 
the $ sector is described by the 0(4) non-linear a model. 
In general, there can be two sources to change the low- 
energy behavior of the non-linear cr model. The first one 
is the existence of a term linear in 9 T $, which follows 
from the analysis, of the equations of motion (EOM) and 
Ward identitiesti As discussed in Ref. [l^, the existence 
of this term relies on the nonvanishing expectation values 



of the corresponding conserved charges, i.e. (J a ) ^ in 
the present case. But (J a ) implies the long range 
ferromagnetic (FM) order and it is forbidden in 1 
dimensions by the Mermin-Wagner-Coleman theorerrc3. 
The second one is the appearance of the 0-term, which 
results from the topological consideration and has no ef- 
fects on the EOM. In our case, £$ does not contain such 
a term because the homotopy group ^[0(4)] = 0. As 
a result, in the long wavelength limit, starts from 
the second derivatives of space and time as shown in 
Eq. (|l9|). Therefore, the broken 0(4) symmetry is re- 
stored and the spectrum corresponding to the $ sector 
is organized as the 0(4) multiplets. In addition, the cor- 
responding excitations acquire an energy gap which is 
given by 



A., = cA, 



(20) 



where c = 0(1) is a non- universal constant. Because the 
charge excitations in this phase are associated with the 
fermion fields and A s <C A e , the mean-field gap A e can 
be identified as the charge gap and the spin gap is deter- 
mined by A s . Although the charge gap is not equivalent 
to the spin gap, there is no spin-charge separation in this 
phase. In conclusion, the long rage excitonic order be- 
comes a short-ranged one and all excitations are gapped. 

A residual Coulomb interaction like the term in Eq. 
(||) leads to the breaking of the 0(4) symmetry and the 
mean-field theory favors & t 7^ 0. Following the simi- 
lar procedure, the Lagrangian which describes the cor- 
responding low-lying fluctuations now consists of Co as 
shown in Eq. (|l8|) and the 0(3) non-linear a model. We 
have no arguments to exclude the possibility of the ap- 
pearance of the 9-tevm in the 0(3) non-linear a model as 
the case of the spin-i Heisenberg chain. If a 0-term ex- 
ists, then the ground state will exhibit an algebraic long 
range triplet excitonic (SDW) order and the spin gap will 
vanish. On the other hand, if the #-term vanishes as the 
two-leg spin ladder, then the spin excitations still have a 
finite gap. 



VI. CONCLUSIONS AND DISCUSSIONS 

In the present paper, we study the low energy physics 
of the undoped nanographite ribbons with zigzag edges 
by neglecting the electron-phonon interactions. We 
show that the interactions between electrons substan- 
tially change the physics of the partly flat bands. Accord- 
ing to the above analysis, there are two possible phases: a 
metallic phase with a spin gap (spin liquid) and enhanced 
SS correlations (region II in Fig. ||), and an insulating 
phase with a gapped spectrum (region I in Fig. [}]). Our 
results that there are two possible phases are valid for 
both the cases of the single layer and AB-stackingtLj. In 
the latter case, the values of v, g, g, and the exponent 
m in Eq. ([ll]) will be different from those in the sin- 
gle layer. Through these parameters v, g, g, and m, the 
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material properties of the graphite ribbon determine in 
which phase this system is truly located. 

To determine which one, the EI, the spin liquid, or the 
gapless edge states predicted by the band theory, is the 
ground state of the graphite ribbon with zigzag edges, 
we suggest two types of experiments: (i) the measure- 
ment of the single-particle DOS by the STM, and (ii) the 
measurement of the uniform magnetic susceptibility. For 
the EI and spin liquid, the single-particle DOS vanishes 
when the energy is smaller than the spectral gap. This is 
distinguished from the gapless edge states predicted by 
the band theory, where the single-particle DOS still has 
a finite value at low energy. Another distinction between 
the SL and the gapless edge states is the temperature de- 
pendence of the uniform magnetic susceptibility. At low 
temperature, it will exhibit the activated behavior in the 
SL due to the spin gap, whereas the Curie-like behavior 
is expected for the gapless edge statescl. 

The enhanced SS fluctuations arise from the singlet 
Cooper-pair tunneling between two bands. It is also the 
origin of the spin liquid phase with strong SS correlations 
in many one-dimensional two-band models as emphasized 
in Ref. [l5|. In contrast, the excitonic fluctuations are 
suppressed after taking into account the quantum fluctu- 
ations. This results from a special feature of the broken 
non-Abelian symmetry in 1+1 dimensions, where the cor- 
responding non-linear sigma model becomes IR unstable 



due to the lack of the topological term. 

In the usual one-dimensional two-band electron sys- 
tem, the competing ordeij-of the singlet superconductor 
is the inter-channel CDWL3, which also results from the 
same interaction as the g term in our Hamiltonian (|Tl|). 
The onset of the inter-channel CDW requires the nearly 
equal electron density in the two bands. This is impossi- 
ble in our case because, here, one band is empty and the 
other is completely filled. 

Finally, we discuss the effects of electron or hole dop- 
ing. For finite doping, the system may be described by 
the one-band Luttinger liquid with a very small Fermi ve- 
locity. For lightly doping, the mean-field theory is sup- 
posed to be robust. Therefore, our results are not af- 
fected qualitatively. According to the analysis of Ref. [I?] 
on the doped EI, there may be an inhomogeneous state 
between the EI and the Luttinger liquid upon increasing 
the doping concentration. 



Acknowledgments 

Y.L. Lee would like to thank M.F. Lin for discussions 
and Global fiberoptics, Inc. for financial support. The 
work of Y.-W. Lee was supported by the National Science 
Council of Taiwan under grant NSC91-2112-M-029-012. 



Electronic address: c001vl00@nchc,gov,t-w 
Electronic address: ywlcc@mail.thu.cdu.t^ 
K. Nakada, M. Fujita, G. Dresselhaus, and M.S. Dressel- 
haus, Phys. Rev. B 54, 17954 (1996). 

K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, Phys. 
Rev. B 59, 8271 (1999). 

Y. Miyamoto, K. Nakada, and M. Fujita, Phys. Rev. B 59, 
9858 (1999). 

K. Wakabayashi and M. Sigrist, Phys. Rev. Lett. 84, 3390 
(2000); K. Wakabayashi, Phys. Rev. B 64, 125428 (2001). 
K. Wakabayashi, M. Sigrist, and M. Fujita, J. Phys. Soc. 
Jpn. 67, 2089 (1998). 

For a review of the excitonic phase, see B.I. Halperin and 
T.M. Rice, Solid State Phys. 21, 115 (1968). 
A discussion of the excitonic phase with the heirachy of 
symmetry breaking SU(4) -> SU(2)xSU(2) -» SU(2)xZ 2 
has appeared in the context of a strong interaction limit 
of the excitonic ordering. See L. Balents, Phys. Rev. B 62, 
2346 (2000). 

For the use of the e-expansion to perform the RG analysis 
in the context of critical phenomena and quantum field 
theories, see, for example, J. Zinn- Justin, Quantum Field 
Theory and Critical Phenomena, Oxford University Press, 
Oxford (1993). 



Maksimov, Sov. Phys. JETP 21, 



9 A.N. Kozlov and L.A. 
790 (1965). 

10 B.A. Volkov, Yu.V. Kopaev, and A.I. Rusinov, Sov. Phys. 
JETP 41, 952 (1976). 

11 L. Balents and M.P.A. Fisher, Phys. Rev. B 53, 12133 
(1996). 

12 H. Leutwyler, Phys. Rev. D 49, 3033 (1994). 

13 In 2 + 1 or 3 + 1 dimensions, (J a ) 7^ is allowed. In that 
case, the term linear in 9 T $ can exist as the FM case 
and the dispersion relation of the Goldstone modes at long 
wavelength will behave as: ui ~ |fc| 2 . 

14 A possible excitonic instability in a related system — lay- 
ered graphite was suggested in D.V. Khveshchenko, Phys. 
Rev. Lett. 87, 206401, 246802 (2001). 

15 V.J. Emery, SA. Kivelson, and O. Zachar, Phys. Rev. B 
59, 15641 (1999). 

16 For a review about the various instabilities in the one- 
dimensional two-band electron systems, see O.A. Starykh, 
D.L. Maslov, W. Hausler, and L.I. Glazman, cond- 
mat/9911286. 

17 L. Balents and CM. Varma, Phys. Rev. Lett. 84, 1264 
(2000); V. Barzykin and L.P. Gor'kov, ibid. 84, 2207 
(2000). 



